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Abstract 

If we consider the gravitational collapse of a material object to a black 
hole, we would expect, for ranges of mass where a black hole would form, 
the following scenario. A large enough object would collapse classically 
until an event horizon forms, and to an external observer the object would 
be lost from view. However, once the horizon has formed, the black hole 
will begin to emit Hawking radiation and the hole will lose mass and 
the horizon will shrink. The final state of this process could be either a 
zero-mass "black hole" with consequent information loss, or some sort of 
"quantum remnant." 

A complete investigation of this process would require: 1) A complete and 
consistent theory of quantum gravity coupled to some kind of field that 
would provide the Hawking radiation 1 (which could be the gravitational 
field itself — gravitons); 2) Some kind of definition of a "horizon" in this 
quantum gravity, and; 3) The calculational tools to achieve a description 
of the scenario. Lacking these, one may resort to toy models to try to 
give some sort of a preliminary answer. 

In this paper we will consider the collapse of an infinitcsimally thin dust 
shell in 2 + 1 gravity, where an exact minisuperspace quantum solution 
exists, and try to make rough estimates of the collapse-Hawking radiation- 
remnant formation process. 



1 Introduction 

There has been some interest over the years in the minisuperspace quantization 
of thin shells as models of the full quantum collapse of more complicated objects. 
One important question that could be asked is whether enough information 
could be squeezed out of the grossly oversimplified model to mock up a first 
look at the complete collapse scenario. In Ref. [2] we used previous work by 
Peleg and Steif |3], Israel [4], and Crisostomo and Olea [5] on the classical 
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equations of motion of a dust shell to construct a quantum formulation of the 
collapse of such a shell. 

This quantum formulation was used to study the possibility of the formation 
of a horizon, that is, whether some sort of "quantum black hole," will form 
or whether a shell of non-interacting particles will simply collapse to a point 
where the uncertainty principle will provide a "repulsive force" and the shell will 
reexpand into the same universe. Simple considerations can answer this question 
up to a point. It seems obvious that for large enough masses one might expect 
a reasonably peaked wave function centered on a radius near to but outside the 
classical Schwarzschild radius of the shell would, as it moves toward zero radius, 
maintain enough coherence to pass beyond the Schwarzschild radius almost in 
its entirety and a "black hole" would form with high probability. For small 
masses one might expect that rapid spreading would overwhelm the coherence 
of the wave function and the shell would reexpand into our own universe with 
probability essentially (or exactly) one. 

In full relativistic quantum gravity these simple ideas are fraught with diffi- 
culties. The most serious of these are: 

1) The quantum shell problem becomes unphysical for masses much above 
the Planck mass. While the minisuperspace approach has frozen out all ra- 
diative modes of the gravitational field, and one is insisting on a single-particle 
interpretation of the shell problem, the wave function in the Schrodinger picture 
can still become pathological at a point where one might expect graviton pro- 
duction to begin. Hajfcek [6] has given a sufficient but not necessary condition 
that shows that we might expect problems above a couple of Planck masses. 
In [7], a qualitative argument was presented that used Compton wavelength 
considerations to give a similar bound. 

2) There are very serious technical problems in the formulation of the prob- 
lem. Many authors have used the full ADM method to construct a Hamiltonian 
for the system, where they have chosen an internal time for the system in order 
to have a true Hamiltonian and a Schrodinger equation that allows the study of 
the time evolution of the quantum system. The choice of an internal time leads 
to the problem of quantum formulations that are not unitarily equivalent. An- 
other problem with these Hamiltonian formulations is that they often require an 
ad hoc choice of a Hamiltonian in terms of variables defined on the shell itself. 
Several of these Hamiltonians have been given by various authors [6], [8]. A 
Hamiltonian due to Hajfcek and Kuchaf [9] has the advantage of being defined 
by a coherent procedure with no ad hoc choices, and is formulated in terms of 
foliations of spacetime by timelike surfaces. We will discuss this Hamiltonian 
in more detail below. All of these Hamiltonians have limits on the mass of the 
shell of a few Planck masses. One formulation that does not seem to have a 
mass limit is based on the Wheeler-DeWitt equation corresponding to one of 
the Hamiltonians in [5]. 

3) Many of the Hamiltonians are quite complicated and there is no real 
chance of finding analytic solutions to the Schrodinger equations of these mod- 
els. Numerical solutions have been presented by a group consisting of A. Corichi, 
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G. Cruz, A. Minzoni. M. Rosenbaum and M. Ryan of the UNAM, N. Smyth of 
the University of Edinburgh, and T. Vukasinac of the University of Michoacan 

The idea of Ref. [5] was to study the quantum collapse problem in 2 + 
1 gravity, where one can address some of the difficulties mentioned above in 
a context where some of the problems mentioned above do not exist. The 
quantum problem is unambiguous for all masses, so there is no problem of wave 
function pathologies. As will be shown below, one possible Hamiltonian for this 
problem has the form of that of a harmonic oscillator. This will allow us to 
find simple analytic solutions that can be used to illustrate the development of 
the quantum collapse of the shell, and allowed us to investigate the problem of 
horizon formation. 

The usual view of black hole evaporation and remnant formation has a vast 
literature. A recent review article [10] cites a large number of authors. This 
usual view has three phases [llj 

1) The balding phase, where the black hole radiates away all its multipole 
moments, losing mass due to classical gravitational radiation. 

2) The evaporation phase, where Hawking radiation of thermally distributed 
quanta carries away mass until the Planck mass is reached. 

3) The Planck phase, where quantum gravity is important and a remnant is 
formed or all the mass is carried away. 

We will not be interested in the balding phase, only in the circularly sym- 
metric evaporation and Planck phases. There are a number of arguments for 
remnant formation (see [10] and references therein), one of which is similar to 
what we will consider here, arguments based on the uncertainty principle. 

Our idea is to investigate Hawking radiation in the quantum collapse of a 
dust shell. This is a minisuperspace model of quantum gravity associated with 
collapse. Our scenario is slight different from the usual view given above, where 
the existence of a stable black hole created some time in the past evaporates 
toward a stable remnant. The back reaction of the Hawking radiation causes 
the mass of the collapsing shell to change, affecting the shell evolution. We 
will attempt to show that this effect causes a shell that has collapsed below its 
(classical) horizon will reappear with a reduced mass. 

One possible way to calculate this scenario would be to define a "shell cloud" 
from the probability density of shell radius, in the same way that we define an 
"electron cloud" from the electron probability density around a nucleus. For the 
electron cloud we can calculate its classical electromagnetic field, while for the 
shell cloud we can calculate its classical metric. The study of Hawking radiation 
in this background could be illuminating. 

Since this idea is difficult to carry out, we will use a rough approximation 
consisting of a classical shell that follows the path of the expectation value of 
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the quantum shell, and a model Hawking radiation. Since in Ref. [2] we saw 
that a shell wave packet collapses to a minimum radius and the reexpands due 
to uncertainty principle effects, this scenario can have the shell falling below its 
classical horizon, and the resulting black hole begins to evaporate, and the shell, 
reappearing with mass loss as some sort of " dynamic remnant" in the form of 
an expanding shell. 

The plan of the rest of the article is as follows. Section 2 will be a brief 
resume of the literature on the general relativistic minisuperspace problem. 
Section 3 will consider the classical and quantum 2 + 1 problems and discuss a 
model of horizon formation, Section 4 will model the complete collapse-Hawking 
radiation-final state scenario, while Section 5 will be conclusions and suggestions 
for future study. 

2 The collapse of thin shells in relativistic quan- 
tum gravity 

The study of the quantum collapse of dust shells is about a decade old [15]. 
The classical collapse problem is fairly straightforward, and can in principle be 
solved exactly. One assumes a J-function massive (or null) shell where Birkhoff 's 
theorem tells us that outside the shell the metric is Schwarzschild and the metric 
inside the shell is Minkowski. The Israel junction conditions can be used to 
derive the equation for the evolution of the shell in terms of intrinsic variables 
on the shell itself, the proper time, t, of an observer riding on the shell and the 
curvature radius, R(r), of the shell that he would measure. The equation for 
the motion of the shell becomes [I] 



where m is the rest mass of the shell (a constant of motion) and M is the 
Schwarzschild mass of the exterior metric. It is straightforward to define x — 
R/m, t — > r/m, V = dx/dr and M/m as our Hamiltonian. Assuming V — 
V(P) and solving dH/dP = V(P) for the "momentum" P, we find that V = 
sinh _1 (P) and our Hamiltonian becomes 



a Hamiltonian given by Hajfcek [6]. Unfortunately, this is not the only Hamilto- 
nian that gives the equation of motion |T]) , and we are left with the problem of 
defining an "appropriate" Hamiltonian for the problem. A number of Hamilto- 
nians for different choices time, including the time of an observer at the center 
of the shell where space is flat, and a Wheeler-DeWitt equation identical to 
that for a relativistic charged particle radially falling in a Coulomb potential 
are given by Hajfcek, Kay and Kuchaf [S]. 




(1) 



H = coshP , 

2x 



(2) 
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Kuchaf and Hajfcek [pj, dissatisfied with such ad hoc Hamiltonians, have 
managed to construct a Hamiltonian for collapsing dust shells that comes di- 
rectly from an ADM reduction of the Hilbert-plus-matter action. The problem 
with this approach is that, as spacetime quantities, the matter variables are pro- 
portional to S[R — Rq(t)], where Rq(t) is the position of the shell as a function 
of shell proper time. It is virtually impossible to reduce the time derivatives of 
these delta functions to reasonable variables in the matter Lagrangian that can 
give a shell Hamiltonian that describes the motion purely in terms of canonical 
variables on the shell. Kuchaf and Hajfcek used an ingenious method (follow- 
ing [T^]) based on the fact that the ADM reduction by a 3 + 1 foliation is not 
restricted to foliation by spacelike surfaces, but works just as well for foliations 
by timelike surfaces. Unfortunately, this approach usually leads to an ill-posed 
problem, but in the case of the shell it does not. Using this approach and a for- 
mulation of the dust fluid velocity in terms of velocity potentials, they define a 
new Hamiltonian. The cost of this consistent formulation is a very complicated 
Hamiltonian, 



Almost all of the Hamiltonians that have been given (except for the Wheeler- 
DeWitt equation of [5]) seem to have mass limits beyond which the wave func- 
tions become pathological. These Hamiltonians are all so complicated that it 
seems impossible to find analytic solutions to assist in the interpretation. In [7] , 
Corichi et al. have given a series of numerical solutions that give the evolution of 
wave functions that are initially sharply peaked over a radius near the classical 
horizon, R — 2M, for the coshP Hamiltonian ([2]) and an approximation to the 
Hamiltonian given by Eqs. (0I3}. All of these solutions show evolution of the 
peak toward R = with a bounce caused by the boundary conditions at R = 0, 
with the appearance of interference fringes as well as a rapid spread of the wave 
packet. The Kuchaf-Hajfcek Hamiltonian has wave functions similar to those of 
(2}, but with many rapid oscillations superimposed. 

Since these Hamiltonians are self-adjoint, unitary evolution implies that a 
peak formed from scattering states will always rebound to R = oo. One can 
ask whether this behavior means that all quantum collapse of this sort implies 
a rebound into our own universe. Since r is proper time on the shell and R is 
also a shell variable, such questions can only be answered by knowing the global 
quantum spacetime surrounding the shell. In any case, the scenario of the shell 
observer is that he sees (begging questions of quantum measurement and the 
reduction of the wave packet) the shell collapse to some point near R = 0, 
where uncertainty principle effects change the classical equations of motion and 
the shell rebounds (actually, a shell where the particles do not interact directly 




(3) 




< R < 2M. 
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with one another "passes through itself" and reexpands, that is, each radially 
infalling particle passes through R = and the azimuthal angle 9 jumps from 
the initial 6q to 9$ + it). Even if the shell has collapsed below its classical 
horizon, in finite proper time it will again be above the horizon and traveling 
toward R = oo. This quasi-classical scenario is not surprising. In proper time, 
a classical shell that manages to avoid forming a curvature singularity at R = 
would behave in this way, but as R becomes greater than 2M the shell would 
be in a universe beyond our temporal infinity (£+), or in "another universe." 

This quasi-classical scenario is what one might expect to see for a large mass 
where the quantum fluctuations would be small compared to the mass, so the 
evolution of the wave packet would be coherent long enough for the shell to 
collapse past its horizon and the shell would emerge from the horizon into a 
new universe. However, if the wave packet spreads sufficiently so that the width 
is greater than the classical horizon radius, we can see that a horizon might 
never form, and the shell would reexpand into our own universe. 

The problem of horizon formation in the quantum system is very difficult. 
Event horizons are global features and one has to try to define a global feature in 
a fluctuating manifold. Of course, this quantum manifold must be constructed 
in terms of the full minisuperspace canonical quantum gravity of the shell- 
metric system. In the shell case we tend to use some kind of approximation 
to construct the spacetime metric. Kuchaf [T^] argues that for the simple shell 
minisuperspace we may just replace the the shell mass (Schwarzschild mass) and 
the shell radius in the metric outside the shell by the corresponding operators 
to make a "metric operator." The problem with this metric operator is that it 
is a function of shell proper time, and studies of the metric close to the shell [17] 
cannot tell us whether a true event horizon (tied to observer time at infinity) 
forms. 

Other approximations are under study [14] . The simplest calculation is to 

calculate < R(r) > and the uncertainty AR. = \J < (R— < R >) 2 > and check 

whether AR becomes very large as < R > becomes small so that < R > does 
not fall below the classical horizon and AR is larger than the classical horizon, 
which can be taken as an indication of the non- formation of a horizon. In [14] , 
numerical evaluations of these two quantities will be presented for the Hamilto- 
nians @ and ([3J3]), and they suggest that no horizon forms for small masses. 
Another possibility (to be considered in [2]) would be to take M\tp(R, r)| 2 to 
be a classical density p(R, r) and calculate the classical metric due to a clas- 
sical fluid with this mass distribution and see whether a horizon forms. Note 
that M should be either the rest mass or the Schwarzschild mass. It is not 
yet clear which. There are technical problems with this calculation. We have 
to calculate the metric from a density that is given in terms of a solution of 
the Schrodinger equation for our Hamiltonian, and there is no guarantee that 
this density can be made to obey the equation T 7 ^ = for our fluid. Another 
possibility considered in [14] is that of the "metric operator" mentioned above, 
which was extended to the whole manifold outside the shell and used to define 
a "quantum" stress-energy tensor. 
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The rotationlcss 2 + 1 problem has some advantages over the 3 + 1 problem. 
The Hamiltonian can be constructed fairly easily, and, as will be shown, has the 
form of a harmonic oscillator. The Schrodinger equation for this Hamiltonian 
has well-known analytic solutions. The expectation value of R and its uncer- 
tainty can, in principle, be calculated analytically. The analogue of the other 
calculation using p(R, r) is much simpler than in the 3 + 1 case. The classical 
horizon is easily found. In the next section these ideas will be considered. 

3 The classical and quantum collapse of shells 
in 2 + 1 gravity 

The first element we need for this problem is an equation for the radius of the 
shell. This problem has been studied in detail by Peleg and Steif [3J , using the 
2 + 1 version of the original formulation of Israel [I] , and Crisostomo and Olca 
[5], using canonical methods. 

Israel studied the collapse of a shell in 3 + 1 gravity, represented by a delta- 
function sphere of dust of radius R(t). In 2 + 1 gravity the shell is a circle, i.e. 
a ring of matter, also of radius R(t). The metric of spacetime will be written 
in circular coordinates, where flat space is represented by the metric 

ds 2 = -dt 2 + dr 2 + r 2 d6 2 . (5) 

The equation for a circle R(t) is 

(3) r = i?(r), (3) 6» = 6», (3) i = t(r). (6) 

We will use the notation i, j = 1,2,3 and A, B = 1,2. We will now need a set 
of coordinates £4 on the circle, which we will take to be £a = (r, 0). 

Off the ring of matter, the 2+1 version of Birkhoff 's theorem says that the 
three-dimensional metric is a static or stationary solution to 

R%j ~ 2 g ij R = -Agy, (7) 

(where we have to have a cosmological constant A to avoid locally completely 
flat solutions) both inside and outside the circle. Here, as in Ref. [2], we will 
study the static case, where the matter has no angular momentum. In this case 
the metric has the form 

ds 2 = -f(r)dt 2 + -^dr 2 + r 2 dQ 2 , (8) 

fir) 

and the well-known solution is 

f = B-Ar 2 , (9) 

B a constant. Since A has units of inverse length, we will write, as is common, 
A = ±l/£ 2 . The final form of the static, circularly-symmetric metric is 

ds 2 = - (-M + ^ dt 2 + j-dr 2 + r 2 d6 2 , (10) 
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where we have taken B, following Bahados, Teitelboim and Zanclli (BTZ) [18 , 
to be —M, M a Schwarzschild mass. We will take the metric inside the ring 
to be the 2 + 1 AdS metric where M = —1 [T5]. [2"0]. As in the 3 + 1 case, 
we expect that outside the ring we will have a black hole metric with some 
"Schwarzschild mass" M, that is, 

ds 2 = - (-M + ^ dt 2 + ^dr 2 + r 2 d6 2 , (11) 

and inside the ring, 

ds 2 = -(l+ r —\ dt 2 + -. 1 2 dr 2 + r 2 dd 2 . (12) 
V PJ (l + £) 

With these preliminaries, in Ref. [2] the equation of motion of the shell was 
calculated using the 2 + 1 analogue of the Israel formulation of the equation of 
motion of a shell in 3 + 1 gravity, which was similar to that used by Peleg and 
Steif[3]. 

In [2] we used the jump in the Einstein equations between the interior and 
exterior of the shell, and the Lanczos relation, 

Iab - 9ab1 = SttSab, (13) 

where jab is the jump of K~ab, the extrinsic curvature of the shell as seen 
from the inside and the outside of the shell, the qab are the components of the 
induced metric on the surface in terms of r and 9, 7 = g AB jAB, and Sab is 
the surface stress-energy tensor, (In our case, we will be interested in a dust 
shell, so we will take Sab = cfuaub, cr the rest mass density of the ring, where 
a = m/2TrR, m the total rest mass of the shell.) to find the classical equation 
of motion for the shell, 

R = -R/e 2 . (14) 
Equation (fT4| has a first integral, 

B-\# + £, (15) 
and, using the Lanczos equation directly, we find 

\/\ + 2E - sJ-M + 2E = 4m, (16) 
which can be solved for E as 

/ M 1 m\ 2 1 

Equation (|15p was used to construct a Hamiltonian formulation of the prob- 
lem. This equation is simply the energy equation for a harmonic oscillator, so 
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it is obvious that if we take E as our Hamiltonian we can define a momentum 
Pr as R, and we have 

Hamilton's equations for this Hamiltonian are equivalent to the equation of 
motion R = -R/l 2 . 

Once we have a Hamiltonian in terms of the ring variables R and t, it is 
possible to construct a Schrodinger equation for the problem. If we write the 
classical relation 

M 1 . . V . - 2 R 2 



— + — + 2m -1 = R 2 + (19) 
8m 8m / t A 

on the right-hand-side M and in are dimensionless, so M = M/Mp\ and m = 
m/Mp\, where Mp\ is the Planck mass, Mp\ = c 2 /G^ 2 ' (we will later need the 
Planck length, Lp\ — hG^ 2 '/c 3 ), where G' 2 ' is the two-dimensional Newton's 
constant with units [L 2 ]/[M][T 2 ]. 

In conventional units, our equation for R becomes 



M Mpi 2m \ 2 1 1 fdR\ 2 R 2 



8m 8m M P , 7 c 2 \ dr £ 2 



(20) 



Our "Hamiltonian" should have the units of energy for a conventional quantum 
Hamiltonian, so we should multiply Eq. (|19p above by Mpic 2 , our i?-momentum 
Pr becomes Mp\R , and with 



M P1 c 2 ( M 2 Mpi 2m \ 2 M P1 c 2 
^-^"18^ + ^ + MpTJ _ ^~' (21) 



and 



"=2^ + i r^ 2 - < 22 > 



where wo = c/i. 

In [5] we quantized this system with H replaced by an operator H. If we 
look at (j2Tj) . we see that the right-hand-side must become an operator, so the 
left-hand-side must also be an operator. We also took M to be a g-number, and 
m and Mp\ c-numbers. Of course, there is no special reason to make this choice, 
but we do so to make contact with previous work (see, for example, |16jV If 
we consider Eq. ©, this choice is motivated by the fact that the Schrodinger 
equation for |(5J) is similar to the hydrogen atom Schrodinger equation, with m 
playing the role of e 2 , and it is usual in the hydrogen atom to take e 2 as a 
c- number rather than a g-number. 

We can now take the wave function of the system to be tp — t{jm'4'(R,t), 
with iPm(M) an approximate eigenstate of M with eigenvalue Mo. An exact 
eigenstate of M of this type would be 6(M — Mo), but to avoid problems with 
the integral of the square of a delta function we will assume that ipM is an 
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extremely sharply peaked wave function centered on M = Mq. In this case, our 
Schrodinger equation becomes (M 2, 0m ~ MqI^m, and realizing Pr as ~iTid/dR) 



MpiC 2 f Mq M] 



PI 



ih 



dr 



8m 8m 
' 2M P1 aR2 



2m V 
MpTJ 



Mhc 2 



(23) 



One curious fact about this equation is that, if we define ip — exp(— iET/K)tj){R), 
then the energy eigenvalues are 



_ , 1 . . , l.Tic 

E n = (n + -)hujQ = {n+ -) — , 



(24) 



and since E is given by the first line of (f2"3")) , there is a discrete relation between 
the Schwarzschild mass, Mq, and the rest mass, m, 



8m 



Ah 



pi 



2 m 



8m Mpi 



-l = (2n + l)^ = (2n + l)-p. (25) 



We will now return to units where G = c = % = 1, Mo now meaning Mo/Mpi, 
m now meaning m/Mpi, and £ meaning i/L-p\. If we solve for Mo in terms of 
m and £ we find 



Mn = 16m 2 



1 



4m 2 2m 2 



1 



(26) 



For a moment we will return to conventional units, and define the follow- 
ing dimensionless variables. We define y = R/y/£Lj>\ and T = ct/1. The 
Schrodinger equation now becomes 



<9T 



2dy 2 



which has eigensolutions 



N /2^( 7 r) 1 /4 



V/a ^n(l/), 



(27) 



(28) 



-ff„ Hermite polynomials. 

In Ref. we studied the evolution of a wave packet sharply peaked around 
a value of y = yo, a point some distance outside the point where a classical 
horizon would form if the radius of the shell were to fall below Rh — iyMo and 
followed its movement as the packet fell toward R = 0. In the 2 + 1 case it is 
possible to give an exact analytic expression for the wave packet as a coherent 
harmonic-oscillator state. However, even though the eigensolutions are nothing 
but harmonic oscillator wave functions, the radial variable R cannot be negative. 
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In order to keep this from happening we will took the potential to be that of a 
half oscillator with an infinitely hard wall at R = This potential is shown in 
Figure 1 



V(R) 




0.5 1 1.5 2 2.5 3 

R 



Figure 1: The harmonic oscillator potential for our problem with Mp\ijJq/2 = 1. The 
dashed line shows a typical position of the classical horizon radius, rn- 

This meant that ip(0,T) was always zero. This boundary condition can be 
enforced by only expanding in odd-n harmonic oscillator eigenfunctions. 

In [2] we began with a difference of two Gaussian states, one peaked around 
y = +yo and the other around y — —yoi so that their sum at y = was zero. 
This state is 

ip(y, 0) = a[e-^ y - yo)2 - e-^ v+y< > )2 }, (29) 

only valid for y > 0, with a a normalization constant. Since this is a sum of 
two Gaussian states, we can use standard techniques to construct the difference 
between two coherent Gaussian states with ([2U]) as an initial condition. The 
result is 

$( y , T ) = ae-^° T/2 e ly o sin 2uJ ° r/i x (30) 
i ex P{-^{(y - Vo cosw r) 2 }) exp(-iyy Q sinc^T") 
~ ex P(-^{(y + Vo coswqt) 2 }) exp(iyy Q sina^r)], 
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which is zero at y = for all r. The normalization a is easily found to be 
a 2 = [VV^(l-e^)]- 1 . 

We connect the variables that describe the wave function with the radius 
of the classical horizon, rjj, shown in Fig. 1. The unnormalized probability 
density, p = vp*ip/a 2 is 

(0 = e -(»-AcosT) 2 +e -( t0 +AcosT) 2 „2e- to2 e- A2cos2T cos(2AwsinT), (31) 

(where w = y/ y/iy/M) . Figures showing the evolution of this probability density 
are given in [2]. This evolution begins with a Gaussian packet at T = that 
collapses toward y — 0, developing interference fringes as it reaches a minimum 
for some y < yo, then rebounds toward y = yo again. This pattern then repeats 
forever. We will only be interested in one cycle of this pattern. 

3.1 The formation of a horizon 

In previous articles the quantization of the shell collapse was used to study the 
possibility of the formation of a horizon in the quantum collapse. As mentioned 
in Sec. 2, this concept has many difficulties. An event horizon is a global 
construct and it has no local definition. This means that in quantum gravity 
one would have to return to the starting point and try to define what a "quantum 
horizon" might be. Once this definition has been decided upon, one must try to 
find out if some collapse process will result in the formation of such a horizon, 
with the result being a probability of horizon formation. Theories of quantum 
gravity in their present state are far from being able to give us this result, 
so, in shell collapse, some articles [7J, [H], [T3] have tried to give an estimate 
of horizon formation by finding out if a sharply peaked wave packet, during 
its collapse toward R — 0, will fall, in some sense, below the classical horizon 
radius, rjy. In some sense, because the packet will usually spread and basically 
will never lie entirely below R = rjj- One could use the integral of ip*ip from 
R = to R = th, which is a number less than one, as the probability of horizon 
formation. Another possibility would be to use the operator M in the expression 

Rh = t\/ 2M — 1 to define a "horizon operator" that could be used to define 
the probability of horizon formation. In previous work it was decided to use 
< R > as a function of r, a quantity that falls from Rq (the R associated with 
yo) to a minimum and then begin to increase again. If this minimum is below 
the classical horizon, one can say that a horizon forms and if not, not. This is 
a yes-no answer instead of a probability, but it is a quick estimate. We will use 
this concept below. 

It is not difficult to calculate < R > (t) from the wave packet given in ([3"T]) , 
but the result is a complicated function that contains the error function and 
Dawson's function, so trying to find the minimum of < R > (r) by finding the 
point where d < R > / dr = requires the solution of a transcendental algebraic 
equation, making it difficult to give an analytic expression for the point where 
a horizon would form. To avoid this problem, in Ref. [2] we used the fact that 
at T = 7r/2, the point where the packet begins to rebound, the peak nearest 
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R = is high and narrow. We used the position of this peak as a parameter 
to tell us whether a horizon would form or not. If the position of the peak is 
below thi a horizon forms, and if not, not. In previous work it was found that 
for large masses the shell collapse was so rapid that the wave packet fell below 
Th so quickly that quantum mechanics did not allow it to rebound before that 
point, while for small masses the rebound occurred for R > rjj. We were able 
to give a range of masses where no horizon would form. 

In the present article, we will study the behavior of < R > (r), but we will 
use the threshold mass where a horizon will form from [2] in order to estimate 
the point where a horizon forms. We will not give any details of the calculation 
in [2], only the results. 

We can define several purely numerical quantities that will be used to explain 
the mass limits. We will key off Th = ly/M, the value of R of the shell as it 
passes its classical horizon. We will also make use of the variable y. The peak 
of the initial wave function is at yo, and using the fact that the value of y 
corresponding to the classical horizon is yn — y/Iy/M, we can define yo — \yu, 
A > 1. The point where the largest peak reaches its minimum at y c — ^yu- 

Studying the motion of the wave packet, it was possible to find an analytic 
relation between A and 7 in terms of Mq, where, for a given A, Mq is a monoton- 
ically decreasing function of 7 (and £). Since 7 must be less than A for collapse 
to make sense, we have a minimum of Mq as a function of £. Another minimum 
can be found by studying H as a function of Mq, m and I. If we take < 2H > 



equal to(f^ + |^- + -p^J - 1, and calculate < 2H > as a function of M , m 
and £, we find a real solution for Mq only for m greater than a minimum value. 
Taking m to be this minimum, we find another minimum for Mq. Equating 
these two minima, we find a numerical value for £, £ 0.18. For this value of £, 
we find that for Mq < 1.48 there is no real solution for Mq, and for 



no horizon forms, while for Mq > 3.32 one does. 

Of course, as has been mentioned above, the rebound of the wave packet, 
once it has passed the horizon does not mean that the shell returns through the 
horizon into the same spacetime where it began. In spite of the fact that in 
terms of proper time the shell exits the horizon in a finite time, in terms of the 
time, t, of an observer at infinity this exit occurs at time after t — 00, or into 
"another universe." 

4 The complete collapse— Hawking radiation— final 
state scenario 

As we mentioned in Sec. 1, we want to consider the complete evolution of a 
collapsing quantum shell. Without a complete, consistent quantum theory of 
gravity, we have to appeal to some sort of approximation. One possible approxi- 
mation that might be the closest we can come to a real theory would be to define 



1.48 < Mq < 3.32, 



(32) 
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a "shell cloud" in the same way we define an "electron cloud" when we consider 
such concepts as electron shielding of the charge of a nucleus. If we consider a 
hydrogen atom with the electron in some quantum state, ip(x,t), we can define 
an electric charge density by p = e4>*ip, and then calculate the electromagnetic 
field due to this "electron cloud." We can do the same with the shell, that is, 
define a classical mass density as something like p = M \ip(R,T)\j i ^ r , and try 
to calculate a classical metric gij (r, r) from this mass density. 

Of course, p(r, t) is not enough to solve the 2 + 1 Einstein equations, we 
need a real stress-energy tensor, T y . If we can define such a T lJ , we can write 
the spherically symmetric metric as 

ds 2 = -e^dt 2 + e^dr 2 + r 2 d9 2 , (33) 

which, if T ij = T^[r, t), we have(see Synge [21]) 

e- x = -Ar 2 + 16vr [ rT$dr, (34) 
Jo 

v = -\-l6ir [ re A (T °-r;)dr. (35) 
Jo 

It is easy to relate T® to the mass density Maip*ip, where a is a numerical factor 
to take into account the fact that we want e~ A to be —Ar 2 — M for large r, but 
Tl is undefined. We also have to study the rest of the components of Tj for 
possible singular behavior. We can show, as in the 3 + 1 case, that T % i should 
diagonal except for T 0r . This means that we need only find T rr and T . The 
Bianchi identities give 

i r Jin) 



T S = --J o r^ d r, (36) 

n = rn + r&-, (37) 

so, in reality, we only need to define as long as Tq can be calculated from 
from Maij)*^. We can consider a fluid stress energy of the type 

Tj = (p + P )u l u 3 + pS}, u e = 0. (38) 

From u l Ui — —1 = — e v (u a ) 2 + e x (u r ) 2 , and defining e x (u r ) 2 = U 2 , we have 

T° = -(l + U 2 )p-U 2 P , (39) 

T; = U 2 p+(l + U 2 )p, (40) 



and 
and 



T£-T; = -(l + 2U 2 ){p + p) 1 (41) 
Tl - p. (42) 
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We can have T ° - T r r = if p = -p and T ° = -p. If we have p = Maip*i;, 
and since we have taken the Schrodinger equation to have the form of Eq. (f23j) , 
what we call ip is actually y/rtp, we will have (taking a = 1/167T in order to 
make e~ x go to — Ar 2 — Mas r — * oo) 

e - A = -Ar 2 - M [ ip*ip(r, r)dr, (43) 
J o 

and f = —A. 

We can calculate Tq and Tjj , which in this particular case (Tq — T£ = 0), 
have singular points. A better choice of U would make Tj well behaved, but we 
have an almost infinite choice of U , although tying U to the quantum mechanical 
current associated with our Schrodinger equation (|23[) would be best. 

Once we have a reasonable metric (taking into account the caveats we have 
mentioned) we could, in principle calculate the Hawking radiation in this back- 
ground, as well as the effect of its back reaction on the quantum evolution of the 
shell. This is still a massive undertaking, involving difficult numerical analysis. 

Our plan in this article is to attempt a simplified analysis, where we can 
use purely analytic constructs (with the exception of a numerical solution to a 
simple first-order ODE) to model the complete collapse scenario. 

The first element we need is some quantity to represent the shell evolution. 
We will consider the expectation value of R, < R > (t) , to represent a classical 
shell evolving with this R(t). We have 

< R > ( T ) = - A ^^ M [cosTerf (\VIVm cos T) 

- sin Te" A2£M ierf {iX^fi^/M sin T)] . (44) 

Figure 2 shows < 7? > (r) for A = 3, I = 0.18, and M = 5. 

This expectation value begins at < R >~ 1.20, falls to a minimum and rises 
again to the original value. In Fig. 2, the value of the classical horizon r# « 0.4, 
is shown as a dashed line, and < R > (r) dips below this value, and later (after 
a relatively short proper time interval) rises above r# again. Of course, the 
amount of time At measured by an observer far from the black hole between 
the moment when < R > (t) dips below the classical horizon and the moment 
when it reappears is greater than infinity. 

The next element needed for the complete scenario is Hawking radiation. 
While Hawking radiation has been calculated for 2 + 1 gravity [55] [23], in 
the spirit of our simplified calculations, we will use a quick estimate of the 
average energy of particles emitted in Hawking radiation and of Th, the Hawking 
temperature. In the Appendix we calculate this energy in 3 + 1 gravity by 
assuming that pair production occurs in regions of the size of the Compton 
wavelength of the particle produced, A c , and that if a pair is produced in a 
small region between r# and A c , one of the pair may fall into the black hole, 
and the other appears as a real particle with velocity c. 
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A simple Newtonian calculation of the energy at infinity of this particle, assum- 
ing that it is created at r# + A c /2 and moves radially toward r = oo, gives an 
energy that is almost equal to the average particle energy in Hawking radiation. 

In 2 + 1 gravity there are several difficulties. One is that there is no Newto- 
nian limit to 2 + 1 gravity, and another is that there is a cosmological constant. 
Since, in the Appendix, we use Newtonian calculations, we can do the same 
here (using two-dimensional Newtonian gravity) . For the cosmological constant 
we can use the Newtonian cosmological constant introduced by Seeliger and 
Neumann [24] at the end of the nineteenth century. 

The two-dimensional Newtonian equation of motion for the radial motion 
of a particle of mass m with a cosmological constant in the spherical field of a 
point mass M is 

d 2 r G (2) Mm mc 2 . . 

m & = — " -pr, (45) 

where c 2 /£ 2 is the Newtonian cosmological constant, which has a first integral 

„ m ( dr\ 2 /on mc 2 9 

E = - - + G^Mmlnr + — r 2 . (46) 



From Figure 3 we have = i\J G^M/c 2 , and A c = h/r 
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Figure 3: The circles showing the horizon th and the thin region between rn and 
th + A c , A c the Compton wavelength of a particle of mass rn. 



and the energy of a particle moving outward at velocity c at a position rjj+A c /2 



is 



r, mc2 „<2u, , ( n Igwm \ c \ i ( £Vgwm x c \ 2 



For X c «r H , 

r 

E « - 



G (2) Mm 




3 VG(2)M 
4*" 



(48) 



The first term is not zero as in the 3+1 problem, but the last term can be 
taken to be the average energy of Hawking radiation particles. This quantity is 
close to the average energy calculated in [22] [23], E = hVG^ M /2tt£, and for 
our approximation we will take it as the average energy. We can now define the 
Hawking temperature (fee the Boltzmann constant), 



3ft VG( 2 )M 



' H 



I 



(49) 



The effect of back reaction can be calculated (using the two-dimensional 
Stefan-Boltzmann equation) as in the Appendix, and leads to a rate of mass 
evaporation of the black hole. The black hole can be treated as a black body of 
radius rn with a surface length, 2-kth, where 



dE 
~dt 



-C(3) 



ftV ( G^M 
i 2 \ c 2 



(50) 
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and the Schwarzschild mass decreases dM/dt = (dE / dt) / c 2 , and, since C(3) (27/32) 
~ 1, we have 



~ -— ( ° mM \ 2 _ f^iiY ( M V 

~df ~ ~e \ c 2 ) =: \TJ \WJ 

or, in our units, 

dM _ M 2 

~dT ~ "l 2 ' 

this can be solved for M(t) as 



(51) 



(52) 



MM M (°) ^ 

-p-t+ 1 

assuming the Af(i) at t = is M(0). 

With these elements we can consider the complete collapse scenario. We have 
to ask whether the classical shell represented by < R > (r) can fall below the 
classical horizon radius rjj = £y/ M(0) and return through the smaller horizon 
at rjj — £y/M(t) at some later, finite, t-time. Naively, one might think that 
the falling horizon radius will meet the rising shell that has passed through its 
minimum radius and is rebounding, "sooner" than it would have in the static 
case. If we were to consider that we have been using some sort of "Newtonian" 
time, which is unique, this would always happen. However, we are trying to 
model relativity, where t should be observer time at large r which is not the 
same as t, the proper time on the shell. 

Classically, we can find a relation between t and r by using the metric pop 
with 9 = 9 = constant., and r = R(t), 

'dT 2 = -(-M+^-)dt 2 + ^-^dr 2 , (54) 

or 



dt sjR 2 + R 2 /e 2 ~M 

Tt |-m + £| ' ( ' 

This expression has the problem that the denominator is zero whenever R = 
t\J~M, so when the shell crosses the horizon (either descending or ascending), 
dt/ dr becomes singular, the origin of the fact that t — > +oo before the rebound- 
ing shell can exit the horizon. 

We want to use Eq. (l55|) modified by quantum considerations with M = 
M(t), M(t) given by Since R 2 + R 2 /t 2 = 2E, we will take this 2E 

to be (2n + l)/£. We now want to use R(t) R > (r). Unfortunately, 

the denominator of (|55)l still becomes zero at < R > (r) = l\J M (t), and 
the relation between t and r is still singular at this point. Arguing that the 
quantum uncertainty in the shell position makes the exact position < R > (r) 
unacceptable in our calculation we can try to model this uncertainty by using 
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< R 2 > [M(t),r] in our expression for dt/tr. Finally, we will try to solve the 
simple differential equation 



(56) 



dr 



with M(t) given by (|53|) . This equation is a nonlinear, first-order differential 
equation that must be solved numerically. 

Our complete collapse scenario will now be: 



1) Collapse from R = \l^J M(0) to the horizon formation point R ~ l^J M(0), 
with essentially no Hawking radiation, since the horizon is necessary for our cal- 
culation of M(t). 

2) From horizon formation to the point where < R > (r) is again above the 
horizon, we solve (|56|) for t{r) and M[t(r)], taking t = at the point where the 
ring crosses the horizon. 

3) Study < R> (t) to see if, on the rebound, it crosses the horizon in finite 
f-timc. 



If this scenario occurs, we will have to discuss the meaning of the "quantum 
remnant" we have. 




0.23 0.235 0.24 0.245 



Figure 4: For M(0) = 5, the time t at large r as a function of r, the proper time on 
the shell. 
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We will now take M(0) = 5 and n = 3 (n must be odd), a somewhat 
arbitrary choice. Note that for n = 3 the maximum M(0) for ([55]) to make 
sense is « 38.9. For A = 3, ^ = 0.18, we can calculate t(r) numerically. This 
was done using Maple, from To where the shell falls below the horizon, to t\ 
where it would reemerge classically with M(0) = constant = 5. The answer is 
given in Figure 4. The use of < R 2 > makes t a regular function of r in this 
entire region. Figure 5 shows the horizon size, fyM(f). 



0.4— 




0.23 0.235 0.24 0.245 

x 



Figure 5: For M(F0) = 5, the horizon radius, th = l\l 'M(r) as a function of r. 



Figure 6 shows < R > (r) and the horizon size t^/M\t(r)} in this region. We 
can see that < R> (r) would never cross the horizon, and we can interpret this 
to mean that Hawking radiation would change M(t) so quickly that no horizon 
would form. 

If we now take M(0) = 10, we can see from Figure 7 that the form of t(r) 
does not change very much. Figure 8 shows t\J M [£(t)] ; and Figure 9 shows 
both < R > (r) and £y/M[t(T)]. In this case, < R > (r) still does not fall 
below the horizon, but a more detailed study of < R > (r) near the initial value 
of r shows that the difference between < R > (r) and r#(r) is only of the order 
of 10~ 5 , so the shell grazes r^(r). 
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Figure 6: For M(0) = 5, the horizon radius th (solid line), and the expectation value 
of the shell radius < R > [M{t(r)}] (dashed line) as functions of r. 




0.225 0.25 0.275 0.3 0.325 



Figure 7: For M(0) = 10, the time t at large rasa function of r, the proper time on 
the shell. 
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0.5- 

r H(*) \ 
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Figure 8: For M(0) = 5, the horizon radius m (solid line), and the expectation value 
of the shell radius < R > [M{t(r)}] (dashed line) as functions of r. 
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Figure 9: For Af(0) = 10, the horizon radius th (solid line), and the expectation 
value of the shell radius < R > [M{t(r)}] (dashed line) as functions of r. 



22 



If we now take M(0) = 20, we can see from Fig. 10 that the form of t(r) still 
does not change very much. Figure 11 shows £y/M[t(T)], and Fi gure 12 shows 
both < R > (r) and tyfM[t{T)] from t ps 0.22139. 
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t(T] 



0.075 




0.05 



0.025— 
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Figure 10: For M(0) = 20, the time i at large r as a function of r, the proper time 
on the shell. 
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Figure 11: For M(Q) = 20, the horizon radius, ru = ty/ M(t) as a function of r. 



23 
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Figure 12: For M(0) = 20, the horizon radius th (solid line), and the expectation 
value of the shell radius < R > [M{t(r)}] (dashed line) as functions of r. 

Figure 13 shows ten times the difference between < R > [M{t(r)}] and rj?(r), 
making the dip below ru more obvious. 




Figure 13: For M(0) = 20, ten times < R > (r) - r H (r), showing the dip of the 
expectation value below th and its return above th- 
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We can see that the shell falls below the horizon at t ps to, ~ 0.22139 and 
reemerges at r ps 0.22241, at which point Hawking radiation ceases. Since i is a 
regular function of r in this region, t at reemergence is finite. We can interpret 
this scenario as loosely defining a "dynamic remnant." Many calculations in 
3 + 1 gravity of Hawking radiation back reaction assume an eternal, static or 
stationary black hole, and that any remnant would be some sort of static or 
stationary construct. In our case, the fact that our quantum shell "rebounds," 
seems to mean that our "remnant" would be an expanding quantum shell. 

The important point here is that there is mass loss during the Hawking 
radiation epoch is « 11, or 55%. For M(0) = 30, the behavior is similar, with 
mass losses of ~ 20, or 67% but the loss is not 100%. As M(0) — > oo, the 
mass loss grows, but it should never be 100%. In the 3 + 1 case, since dM/dt 
decreases with growing M, we might expect the mass loss to decrease. 

Our scenario is the collapse below a horizon, mass loss due to Hawking 
radiation, and the reappearance of the shell in expansion, but with a reduced 
mass. We can consider this a "dynamic remnant." There are a number of 
problems with this calculation beyond its rough nature. One is that we have 
assumed in our Hawking radiation calculation that the Compton wavelength, 
A c , is much smaller that r# . However, our units show that M is of the order of a 
few Planck masses, so A c for any reasonable particle is very much larger than r# . 
We would have to drastically modify our Hawking radiation approximation in 
this case. This problem reflects the well-known difficulty of calculating Hawking 
radiation for particle masses above the horizon radius (see [10] and references 
therein). 

Another problem is that we can still have information loss. The existence of 
a remnant does not necessarily solve all problems. There is always the necessity 
for an infinite number of states which allows for the unbounded information 
content inherited from the original state. 

5 Conclusions and suggestions for further re- 
search 

While the above results are suggestive, there are several caveats, and one should 
be careful in interpreting the results. One problem is that Eq. (|52p has the rate 
of mass loss increasing as M increases, while in 3 + 1 gravity it decreases, 
so even if we accept the results, they will be numerically quite different from 
those in 3 + 1 gravity. Of course, our results are, at best, crude, due to the 
drastic simplifications we have made. Perhaps the most sensitive choices we 
have made are in the equation for dt/dr. If we were to use the approximation 
of the classical T y constructed from ip*ip, some of these difficulties could be 
sidestepped, but there will still be a difficult problem in that the proper time r 
is valid only for one shell, and we would somehow have to change the problem 
to account for different proper times on different shells. One advantage we have 
in our simplified calculation is that we can argue that our proper time is valid 
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for the one shell we have. 

Another time problem is that we are using (essentially) Schwarzschild time 
for t, while a coordinate system such as Kruskal might be a better choice. 
However, the coordinate change from Schwarzschild contains M, and a changing 
M makes the coordinate change difficult. In the classical T y problem, we would 
have to invent the equivalent of Kruskal coordinates. 

The one advantage of 2 + 1 gravity is that the quantum collapse problem is 
exactly solvable, but numerical work in 3 + 1 would be more believable. 

However, we can make several suggestions for further research: 

1) Do the T y problem, with due attention to the difficulties with proper time. 

2) Return to 3 + 1 gravity and try to do a similar simplified calculation to 
the one attempted here. 

3) Redo the calculation with due care given to the fact that, in general A c >> 
Each of these would require numerical solutions. 

Appendix 

We want to do the same simplified calculation of Hawking radiation energy that 
we did in Sec. 4 for the 2 + 1 black hole. If we use Fig. 3, taking the circles 
to be spheres, we have a thin shell between th and tr + A c where we can have 
the production of virtual particles with mass m, with one of the pair falling 
into the hole, and the other becoming real with outward radial velocity c. The 
(constant) Newtonian energy of the particle, if produced at r — Rh + A c /2, is 




th + A c /2 



GMm 



(57) 



and if A c << rjj, 




(58) 




(59) 



2 2 8GM 8GM' 
The exact average energy of Hawking radiation particles at infinity is 




(60) 



which is within 14% of our simplified result. 
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We can now define a Hawking temperature 

he 3 
8GMk B ' 



-H 



(61) 



and treating the black hole as a black body with radius Th, we have an energy 
loss due to Hawking radiation as (<r the Stefan-Boltzmann constant) 

§ = -a(Anr%)n, (62) 



15360 \G 2 J \M 2 J 

The mass loss equation is 

dM _ 7T 3 /tu?\ f 1 \ 



(63) 



dt ~ 15360 \ G 2 ) \M 2 )' (64) 

Note that dM/dt decreases with increasing M, rather than growing as in the 2 
+ 1 case. Also, Tu is independent of m, so, formally, we can take m — > to 
assume the production of massless particles. 

Notice that we have assumed Th « A c , and if we wanted to use our sim- 
plified calculation for Planck-sized black holes, we would have to change our 
approximation. 
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